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ABSTRACT 

Although several dozen double white dwarfs (DWDs) have been observed, for many the exact nature 
of the evolutionary channcl(s) by which they form remains uncertain. The canonical explanation calls 
for the progenitor binary system to undergo two subsequent mass-transfer events, both of which are 
unstable and lead to a common envelope (CE). However, it has been shown that if both CE events 
obey the standard acE-prescription (parametrizing energy loss), it is not possible to reproduce all of 
the observed systems. The 7-prescription was proposed as an alternative to this description, instead 
parametrizing the fraction of angular momentum carried away in dynamical-timescale mass loss. In 
this paper, we analyze simultaneous energy and angular-momentum conservation, and show that the 
7-prescription cannot adequately describe a CE event for an arbitrary binary, nor can the first phase 
of mass loss always be understood in general as a dynamical-timescale event. We consider in detail the 
first episode of mass transfer in binary systems with initially low companion masses, with a primary 
mass in the range 1.0-1.3 M Q and an initial mass ratio between the secondary and primary stars of 
0.83-0.92. In these systems, the first episode of dramatic mass loss may be stable, non-conservative 
mass transfer. This strips the donor's envelope and dramatically raises the mass ratio; the considered 
progenitor binary systems can then evolve into DWDs after passing through a single CE during the 
second episode of mass loss. We find that such a mechanism reproduces the properties of the observed 
DWD systems which have an older component with M < 0.46 M Q and mass ratios between the 
younger and older WDs of q > 1. 

Subject headings: binaries: close — stars: evolution — X-rays: binaries 



1. INTRODUCTION 

A critical phase in the formation and evolution of 
most close interacting binaries is the so-called common- 
envelope (CE) event, during which the components of a 
binary system are engulfed by a common gaseous enve- 
lope. The resulting drag-like interaction wil l then dra- 
matic ally shrink their orbit (IQstrikerl 119761 : IPaczvnskl 
Il976h . Depending on the envelope structure and com- 
panion masses, either the envelope is ejected during this 
process leaving behind a close binary, or the two stars 
merge. Dividing the parameter space into binaries that 
survive CEs and those that do not, and determining the 
final separation for the former, is necessary in order to 
calculate the formation rates of low-mass X-ray binaries 
(LMXBs ) and 7-ray bursts, as w ell as for LISA and LIGO 
sources (|Belczvnski et all I2008D . Understanding which 
evolutionary channels are viable pathways for the cre- 
ation of double-degenerate systems such as double white 
dwarfs (DWDs) is one of the key 'sanity checks' avail- 
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able in studying binary evolution, providing tight ob- 
servational constraints for the otherwise almost entirely 
theoretical study of the evolution of multiple stars. 

The majority of currently observe d DWDs consist 
of re l atively low-mass helium d warfs (INapiwotzki et aT] 
2001 [200a IKarlet alJ 12001 INapiwotzki et al.l 120041: 
Nelemans et all 120051 ). The observed DWDs for which 
the older companion is less massive than 0.46 M Q have 
an average orbital period of 0.96 d and an average mass 
for the (inferred) younger companion of 0.39 M Q (based 
on the doubl e-lined spectroscopic bina ries observed prior 
to this year. Ivan der Sluvs et al.ll2006l see table 1). 

In order to form WD companions of the observed 
masses, a DWD system must necessarily have undergone 
two episodes of envelope mass loss as each component 
went through the giant phase. This can be inferred from 
the fact that none of these observed WDs are sufficiently 
massive to have evolved independently to a degenerate 
state within a Hubble time, and so their progenitors, 
though still low-mass stars, had to be several times more 
massive than the WDs they produced. As low-mass gi- 
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ants obey a tight relation between core mass and radius, 
the mass of the younger WD gives the strongest known 
constraint on the pre-mass-transfer binary configuration. 
On the other hand, the initially lower-mass companion 
must evolve into a white dwarf via a common-envelope 
process in the last binary interaction ph ase, in order 
to achieve the small s eparations observed (jlben fe Livid 
119931: lWebbink|[2008l ). However, the physics underly- 
ing the evolutionary channel(s) by which DWDs form 
remains uncertain. 

A variety of explanations have been posited to explain 
the first stage of mass transfer (MT), in w hich the ini- 
tial primary loses its envelope. [Han (1991) showed that 
both CE and Roche-lobe overflow (RLOF) are possible 
in the first episode of the MT. Later, it was claimed that 
this ph ase cannot be described by stable and conservativ e 
RLOF (jNelemans et al.l[2000t Fran der Sluys et alJl2006T >. 
The suggestion that the first MT episode in the formation 
history was a CE event also proved problematic: when 
formulated in terms of the binding energy of the enve- 
lope, it appeared to imply that the binary must be able 
to eject it s envelope with an unphysical (often negative) 
efficiency dNelemans et al.ll200Ct TNelemans fe Toutll2005t 
iWebbinkl |2008|). Furthermore, the majority of the dis- 
covered double-degenerate systems have components of 
comparable mas s, in contra diction to the results of past 
numerical work (|Hanlll998h . 

Recently, a reconciliation was attempted by in- 
stead considering angular-momentum balance, re- 
parametrizing the problem in terms o f the so-called "7- 
prescription" (|Nelemans et al.l 120001 iNelemans fc Tout! 
I2005t ). In this paper, we demonstrate that the 7- 
formalism is not stable against small changes in the bi- 
nary parameters (§3), and that in general it too fails to 
provide a physical description of an initial, dynamical- 
timescale phase of evolution which is consistent with 
both energy and angular momentum conservation. 

As an alternative, we remove the restriction to 
dynamical-timcscale processes and consider the evolu- 
tion of red-giant main-sequence (RG-MS) binaries via 
stable, but non-conservative mass transfer from the pri- 
mary as the first phase of mass loss (§ 4, 5). In this paper, 
we limit ourselves to systems where the initial primary 
has not reached the helium flash. We demonstrate that, 
upon the former secondary reaching RLOF, an ensuing 
CE phase will produce a set of DWD binaries with mass 
ratios and periods in line with the observed DWDs whose 
older remnant is < 0.46 M Q (§6). Therefore the stable 
MT+CE channel provides a natural means for reproduc- 
ing those DWDs in which the first phase of mass loss 
appears to have been accompanied by orbital expansion. 

2. PRELIMINARIES 

We now outline the two canonical prescriptions for 
treating common-envelope evolution, and develop the 
standard formalism for treating stable mass transfer in 
the case of Roche-lobe overflow. 

2.1. Energy balance versus angular-momentum balance 

In the standard treatment of CE outcomes, the final 
separation of the binary is determined via the "energy 
formalism" ( Wcbbink 1984), in which the binding energy 
of the (expelled) envelope is equated to the decrease in 
the orbital energy E ^: 
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Here a; and at arc the initial and final binary separa- 
tions, TOd and m a are the initial star masses (donor and 
accretor, respectively) and md, c is the final mass of the 
donor, after losing its envelope. 

A parameter A is introduced to characterise the central 
concentration of the donor envelope: 
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Here md, G is the mass of the removed (giant's) enve- 
lope, r d is the radius of the giant star at the onset of the 
CE and e is the specific internal energy. -Ebind therefore 
consists of the potential energy of the envelope and its in- 
ternal energy, and can be found directly from the stellar 
structure for any choice of core mass. In our calcula- 
tions, we adopt to include in e only the thermal energy 
of the gas and the radiation energy, but not the recombi- 
nation energy. There are also altern ative definitions for 
-E-bind, where ionizati on energy (e.g. lHan et all I2001 or 
enhanced winds (e.g. Soker 2004) arc taken into account. 

Another parameter, «ce, is introduced as a measure 
of the efficiency with which energy is transferred from 
the orbit into envelope expansion. Invoking energy con- 
servation, one can then find the final orbital separation 
from 
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The obvious bounds on this parameter are < «ce < E 
since energy cannot be generated. The parameter A 
can be calculated from detailed stellar-structure mod- 
els. For many low- mass giant stars of interest A ~ 1, 
and therefore many auth ors assume aX = 1 or 0.5. 
Ivan der Sluys et all pOlOl ) find 0.027 < A < 1.73 with 
a median of A = 0.86 for the donors (with masses up 
to 20 Mq) of a population of 165,000 binaries at the 
onset of a CE. Note that this do es not hold for mas- 
sive giants, where 004 < A < 10 (|Tauris fc Dewil I200H 
IPodsiadlowski et al.|[200l ). 

An alternative approach to determine the orbital pe- 
riod after the ejection of the envelope is by considering 
the other available integral of motion, the angular mq - 
mentum (|Nelemans et alj|200fl INelemans fc Toutll2005t) . 
This is known as the 7-formalism and can be expressed 
as: 
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where Ji and Jf are the angular momenta of the ini- 
tial and the final binaries. Hence, the 7-parametrization 
describes the specific fraction of the initial angular mo- 
mentum that is carried away by the envelope as it is lost 
from the system, in terms of a multiplicative factor (7) 
times the fraction of the total mass lost from the system. 

It was shown that the 7-formalism can reproduce the 
distributions of the orbital periods and mass ratios in the 
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observed DWD systems with a single value of 7 = 1.5 for 
the first mass-transfer event, where the range for possible 
recons tructed values was found to be between ~ 1.2 and 
~ 2.5 (jNelemans fc Toutl[2005h . At the same time, they 
showed that no single value of oiceA can satisfy all the 
merger scenarios they considered. The values that can re- 
construct the observed systems scatter from unphysically 
negative to suspiciously large ~ 2, with values for «ceA 
skewed to low (< 0.5) values. We note that those low 
values are rather in agreement with expected «ce < 1 
and A variation along the giant branch. The considered 
systems had primary masses from 1.2 to 3Mq, and most 
of the likely progenitors clustered around ~ 1.7 Mq. 

2.2. Orbital evolution during the mass transfer 

Assuming a circular orbit, the instantaneous angular 
momentum about the binary's centre of mass is defined 
as 



J = 



(5) 



where a is the current orbital separation. 

Neglecting the spin angular momentum of the com- 
panions, a general equation for the orbital evolution of a 
non-eccentric binary is: 
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Here J and a are the rates of change of the orbital angular 
momentum and of the orbital separation, respectively. 
The quantities m d and m a are the corresponding rates 
of the total mass loss from the donor and the accretor. 

The total mass loss from the donor is rn d = ihd,t + 
Wd,wj where TOd,t is the mass transfer rate through L\, 
and ?fid,w is the mass loss from the donor in the form 
of a fast isotropic wind. As the accretor is less evolved 
than the donor and is still on the main sequence, its wind 
mass loss can be neglected. We introduce the following 
notations: 

• = — m a /md,t — the fraction of the transferred 
mass that is accreted onto m a ; the conservation 
factor. With — 1 the MT is fully conservative. 

• 8 = m^w/md — the fraction of the total donor 
mass that was lost in the form of a fast isotropic 
wind. 

With these notations, the total mass loss from the sys- 
tem is 

m d + m a = (1 - + 08) ma, (7) 



and the mass accreted is 

m a = —0(1 - 6) m d . 



(8) 



The specific angular momenta of the two binary mem- 
bers, and hence of the mass that is carried away from 
them (from the accretor star as transferred material that 
is isotropically re-emitted, and from the donor star in the 
form of a fast isotropic wind), are: 
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and the total loss of angular momentum from the binary 
system can be written as: 
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Combing eq.® and ca. ([TUl) . we find that the orbital 
evolution for our calculations of evolutionary sequences 
during RLOF is: 
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Note that in our calculations we only choose to fix the 
constant conservation factor for each individual binary 
evolution, but not 8, which is always a function of the 
current evolutionary state and MT rate. During the ini- 
tial MT phase, when MT often operates on a timescale 
close to thermal, and the MT rate is high, 8 is effectively 
0. The inclusion of 6 becomes more important during 
the late stage of the MT, when the MT operates on the 
nuclear timescale (for more details, see § 5). 

3. ANALYSIS OF THE 7-FORMALISM 

iNelemans et al.l (|2000D analyzed the formation of 
DWD by attempting to reconstruct known systems. The 
authors found that for the second phase of mass loss, 
a CE event with «ceA ~ 2 could be fit to the ma- 
jority of observed systems. However, one should assess 
the implications of this with caution as we would expect 
use of a single-valued parameter to be inappropriate in 
the general case, as there is little reason to assume that 
such a parameter should not vary under differing condi- 
tions. Specifically, in the case of the acE^ formalism, A is 
known to vary along the giant branch, and it would seem 
a somewhat awkward contrivance to suggest a single en- 
velope removal efficiency is shared among all binaries un- 
dergoing a CE. Such considerations are essential in order 
to prov ide more substantive r esults in population syn- 
thesis. INelemans et al.l (|2000f ) find plausible values for 
ckceA ranging from 0.5 to 3 for different observed cases. 
We take a value of acE to be reasonable so long as it is 
in the physically meaningful range betw een and 1. 

The same reconstruction analysis in INelemans et al.l 
(|2000f ) has also shown that the first phase of mass loss 
can be explained with a CE only if —15 < «ceA < —5, 
well outside any physically acceptable range. However, 
this assumes that the only relevant source of energy is the 
binding energy of the envelope. For the orbit to widen 
in the course of mass loss, Eq.© requires an influx of 
energy outside of that which is available on a dynamical 
timescale. This is not at all unphysical for dynamically 
stable MT via RLOF, which may take place on the ther- 
mal or nuclear timescale and is driven by the nuclear 
burning of the donor. Therefore the cause for this un- 
physical behaviour would seem to be immediately clear; 
stable MT via RLOF can in many cases provide a clear 
solution for the removal of the initial primary's envelope. 

However, the possibility that stable conservative MT 
could account f o r this stage was previously ruled out. 
INelemans et al.l (|2000l ) empirically analyzed Algol type 
MT, specifically MT initiated while the donor is still on 
the m ain sequence. In a later study (|van der Sluvs et al.l 
I2006T ) . conservative MT was considered and rejected once 
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more. We note however that neither of the studies per- 
formed the deta iled evolution fo r gian t donors at the 
start of the MT. iNelemans et al.l (|2000T ) empirically an- 
alyzed only stars with less than 50% of the outer enve- 
lope being convective at RLOF (e.g., for a 1.2 M© star, 
this limit means that the star is as small as ~ 2.6 Rq 
and has barely left t he bo ttom of the giant branch). In 
Ivan der Sluvs et all (|2006l ). it was adopted that all stars 
beyond the bottom of the giant branch will start un- 
stable MT (see their condition for instability 7), so the 
detailed MT calculations were performed only for stars 
that have not become giants yet upon reaching RLOF. 
We note that in both studies the donors could become 
giants after the start of the MT and continue MT stably, 
but the initial conditions limited the donors at the end 
of the MT to rather early giants. 

An alternative formulation, in which the envelope of 
the initial primary is lost without significant spiral-in, 
was formulated in terms of a new parameter 7 (see W2.1\i . 
It was found that 7 = 1.5 allows for widening or very 
mild shrinkage in the first CE phase (for mass ratios very 
close to unity) but extreme sh rinking in the second. As 
a result, it has been claimed in INelemans &: Tout! ()2005l ) 
that this single value of 7 can explain all observations 
and therefore can be used to predict the outcome of a 
CE, as well as any other dynamical-timescale evolution, 
for any arbitrary binary population, including massive 
binaries. 

The 7 formalism can indeed lead to orbital expansion, 
as well as an increase in the total energy, for mass ra- 
tios close to 1. However, if interpreted as a dynamical- 
timescale process, this increase must prove just as un- 
physical as a negative efficiency in the oceA formalism. 
In the following section we will also demonstrate that a 
single- valued 7 is not robust against small perturbations 
in orbital parameters, and is thus unable to explain both 
evolutionary stages in all instances. Therefore, when at- 
tempting to describe both phases of DWD formation si- 
multaneously, the 7 formalism is subject to similar funda- 
mental physical problems as those found using the ckceA 
prescription. Rather, in § [5] we demonstrate that stable 
but non-conservative mass transfer can account for the 
expansion in the initial phase, bypassing the need for any 
parametrization of this stage. 

In order to understand what the 7 mechanism im- 
plies, and to contrast it with the ckceA formalism in 
a double- CE scenario, let us analyze the energy and 
angular-momentum balance during a first and second 
CE of a binary system that is a potential progenitor 
of a DWD system. After the first CE, this system will 
presumably evolve to a second (and last) CE event, ac- 
co rding to the double- CE model via the 7 mecha nism 
bv INelemans et al.l (|200l : INelemans fc Toutl ([20051 ). In 
their scenario, at the start of the second CE event a giant 
donor has a mass of about 1.7 — 2.2 M© and the original 
primary star has already formed a WD with a mass of 
- 0.35 - 0.55 M Q . 

3.1. Common envelope via a-formalism 

As an example, we consider a binary with a giant of 
2M© and a WD of 0.5 M , where the CE is initiated 
when the giant has a radius of 21 Rq. The states of the 
binary system in the orbital-energy-angular-momentum 
(E — J) plane are plotted in Fig.[TJ The location of the 
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Fig. 1. — Orbital angular momentum J and energy E for a CE in 
a 2M Q +0.5M q binary (r d = 23_R and m djX = 0.337 M ). The 
black star indicates the state of the binary at the onset of the CE. 
The solid line (red) shows the Kcplcrian -E or b — J rb relation for the 
final binary assuming that the mass of post-CE remnant consists 
of all the RG mass that was originally contained within the final 
Roche lobe (mm), in other words, the maximum possible remnant 
mass for this orbit. The shaded region contains Keplerian orbits 
for core masses bounded by m c m ; n < m d c < m c ,max (see the 
text). A final state of the post-CE binary must lie within this area. 
The dashed lines (blue) trace the minimum energy expenditures to 
release the envelope with, for A calculated using the cq. {2} and for 
A = 1. Dashed-dotted green line separates regions where the post- 
CE binary has more or less energy than at the start of the CE. 
Dotted green line separates regions where the post-CE is wider 
(a > a;) than at the start of the CE and where it shrunk (a < a;). 

initial state is marked with a star. For this state the or- 
bital energy -Eorb.i and the rotational energy of the giant's 
envelope (assumed to be synchronized with the initial or- 
bit), as well as the orbital angular momentum J or b,i and 
angular momentum of the giant's envelope are taken into 
account. 

The thick solid line is determined by repeatedly re- 
moving the outermost mesh point from the structure 
model, assuming that the radius of the resulting star is 
equal to the radius coordinate the new outer mesh point 
had in the original model (i.e., the model doesn't ex- 
pand), equating this radius to the Roche- lobe radius and 
then computing the orbital properties of the resulting 
binary. We also note that in this figure, the adiabatic 
expansion of the post-CE remnant due to mass loss was 
not taken into accou nt (see discussion in lDelove fc Taaml 
l2010t ll^novd[20Tl| ). as such, the solid line provides the 
closest possible orbits and, realistically, the position of 
the binary can only be below this line. The region where 
the energy is below its initial value is an impossible final 
configuration for a post-CE binary, and is present in the 
figure only as a limiting case, as such a configuration is 
not physically possible. 

For any given core mass, Kepler's law constrains the 
orbits to a straight line parallel to the shaded region, 
which is bounded by the minimum and maximum core 



The formation of double- WD binaries 



■5 



masses. Our choice of the lower bound on the possible 
core mass m Cim ; n is the hydrogen-exhausted core m^x 
(the region where X < 10~ 10 ) and is rather standard. 
For the upper bound on the possible core mass m c ,max 
we choose the minimum between the cen tral mass which 
conta ins less than 10% hydrogen (as in IDewi fc Taurisl 
120001) and the bottom of the outer convective envelope 
("^bce)- The definition of the post-CE core is an im- 
portant problem in itself. We should stress that there is 
no strong constraint from either observations or theory 
which dictates that the post-CE remnant immediately af- 
ter the dynamical event will have the exact same mass as 
the WD companion that is observed much later. As such, 
the mass of a post-CE remnant could be anywhere be- 
tween our minimum and maximum core-mass definitions. 
A stron ger constraint on the po s t-CE c ore mass ^pro- 
posed in llvanova fc Chaichenetsl ()2009f ): llvanoval (|2011[ ). 

The solid line in Fig.[T] represents the closest possible 
orbit from Roche-radius considerations: only the states 
below the solid line are possible, and the actual state 
depends on the degree of the post-CE detachment. Only 
the strip bounded from above by the thick solid line and 
lying within the shaded region is permitted for a self- 
consistent final post-CE binary. We note from the FigQ] 
that, depending on the choice of the post-CE remnant 
mass, the minimum possible semi-major axis of the post- 
CE binary can vary by a factor of 30. 

The post-CE binary can not have more energy than the 
initial binary and still satisfy conservation of energy (see 
dashed-dotted green line on Fig. [1]), since some energy 
must be used to expel the common envelope to infinity. 
We can also find the position that a final binary would 
have if it kept the same orbital separation, for any valid 
post-CE mass of the donor m c>m i n < md, c < n~*-d- Clearly, 
a widened orbit would violate energy conservation. 

We also plot the energy "expense", the initial energy 
value plus the energy needed to disperse the envelope 
above a given core mass, as the dashed lines in Fig.Q] 
(assuming a = 1 — the most efficient case possible, a 
smaller value would shift the curve higher). The differ- 
ent points on these tracks are determined by removing 
the outer mesh point, computing the amount of energy 
needed to expel that mass shell to infinity, and removing 
that energy from the binary orbit. The angular momen- 
tum for these states is computed using the remaining 
masses and orbital energy to satisfy Kepler's law. As 
this energy expense is the smallest possible, the final bi- 
nary state should lie anywhere above the dashed lines 
for a given core mass. For comparison, we show a solu- 
tion where A is calculated from the stellar structure, and 
where A = 1 is assumed. 

Energy considerations limit the range of possible post- 
CE binary separations, as there is a minimum amount 
of energy that needs to be extracted. Nonetheless, for a 
properly computed A, final orbital separations can vary 
over a factor of 10, depending on the choice of core mass. 
The whole range of possible states is now bounded by the 
blue dashed line to the bottom, the solid red line to the 
top and the shaded red area. When the CE is initialised 
at a different giant radius (i.e., at a different orbital pe- 
riod), and accordingly a different m^x, the picture is 
qualitatively similar, although the uncertainty that is in- 
troduced by the core-mass definition could vary. 



3.2. Common envelope via ^-formalism 

We now construct an exa mple that fits t he 7 - 
parametrisation introdu ced in | Nelemans et al.l (1 20001 ) 
and further studied in iNelemans fc Toutr ~( 2005[) . and 
show that it is not robust against small changes in bi- 
nary parameters. The latter study determined the range 
of values for 7 to be between ^1.2 and up to 4 in the 
first episode of the MT and between ^0.5 and 3 in the 
last episode of the MT, with the most likely value of 
7 m 1.5 — i.e., every gram of the escaped envelope ma- 
terial carries away 1.5 times the average specific angular 
momentum of the binary. 

Let us choose a 2M giant, evolved to a hydrogen- 
exhausted core mass m,i.x ~ 0.4M fl , a common case 
in the study of INelemans fc Tout! (|2005[ ) , where many 
double WD binaries have an older companion of 0.5 M Q 
and a younger WD of 0.4 M@. In Fig. [5] we again plot 
the final binary configurations in the E — J plane, in 
accordance with the 7-formalism. The black tracks arc 
constructed by subsequently removing mass-mesh points 
from the surface of a stellar-structure model. The value 
for J is computed from the original angular momentum of 
the binary, and by applying the 7-formalism. The value 
for E is computed then from Kepler's law. The only self- 
consistent solutions are those that fall within the strip 
delineated by the two thin solid red lines (corresponding 
to the red shaded area in Fig.Q] and described in § 13. 1|) 
for the adopted core-mass range. The solutions also are 
bound from above by the solid red line, where the binary 
is bound to merge, and by the dashed-dotted green line, 
which denotes the initial binary energy: below this line, 
the 7-formalism essentially predicts energy generation. 

From Fig. [21 we see that the track for 7 = 1.5 roughly 
coincides with a possible final binary configuration for 
these particular companion masses. By coincidence, the 
7 « 1,5 solution crosses the final binary configuration 
at approximately the location mandated by energy con- 
sideration. However, we will show below that this is no 
longer the case for a giant of the same mass but with a 
different core mass. The self-consistent solutions in this 
example are those with 7 « 1.43 — 1.47 (dotted black 
lines), where for some core masses values between 1.38 
and 1.505 could also give physically valid solutions (thick 
solid black lines). Note that the change in 7 by ±0.05 
around 7 = 1.45, for the same core mass, results in final 
orbits being different by a factor of 10, while a bigger 
change in 7 leads to either merger, or to an unphysical 
energy generation. As such, it is not possible for one 
system to have a large range of 7 givi ng physically valid 
solutio ns, such as from 1.2 to 3, as in INelemans &: Tout! 
(|2005T ). 

As we indicated, the consistency of the solution con- 
sidered above is coincidental. We will now consider two 
examples where a CE with 7 a 1.5 does not lead to a 
self-consistent solution. 

First, we take the same 2M Q giant with the same 
TOc.min ~ 0.4 Mq, but with a 1.5 M companion (Fig.[3|), 
similar to some ini tial binaries considered to be progen- 
itors for DWDs in Ivan der Sluvs et al.l (|2006h . and con- 
sider how the 7 formalism treats the first phase of MT 
here. This binary would have a self-consistent post-CE 
binary separation only for 7 w 1.99 — 2.06 (close to values 
of 7 used in some similar systems in Ivan der Sluvs et al.l 
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Fig. 2. — Orbital angular momentum and energy for a 2Mq + 
0.5 Mq binary (r d = 86i? and m d , x = 0.38 Mq). Thick and 
thin solid red lines show the only possible final binary configura- 
tions for various adopted core masses, for the comparison we also 
show the the minimum energy expenditures to release the envelope 
(blue dashed line, see also the caption of Fig.[TJ. Black solid and 
dotted lines indicate possible final binary configurations, assuming 
angular momentum lost in accordance with the 7-formalism, where 
thick black solid lines show minimum and maximum 7 that make 
a physically viable binary, and dotted black lines show values of 7 
that lead to formation of physical viable binaries for all range of 
the core masses. Other lines' styles are as in Fig.[l] 

(|2006ft . see their Table 6). With 7 = 1.5, this binary 
would extract too little angular momentum to remain an 
energy-consistent binary: the binary is even less bound 
than it was at the beginning — the same effect as hav- 
ing negative a or a stellar wind; apparently energy is 
generated during the mass removal to infinity. As a con- 
sequence, the binary is becomes wider during the mass 
loss. Speaking physically, it means that, if 7 = 1.5 can 
reconstruct some systems, in these systems the mass loss 
process was not conserving energy, and as such could not 
proceed on a dynamical-timescale and be called a com- 
mon envelope. This is a clear indication of a stable mass 
transfer event and the MT phase has to be studied and 
understood using the appropriate tools. 

Another example involves a binary with a smaller com- 
panion (0.35 Mq, Fig-HJ) - Kcplcrian solutions can be 
found for 7 w 1.36 — 1.39, where the largest acceptable 
solution is 7 ps 1.415. With 7 = 1.5, this binary has 
to extract so much angular momentum that the post-CE 
binary clearly has to merge if the realistic core sizes are 
taken into account, whereas energy consideration indi- 
cates that a merger is not necessary, though it is at a 
threshold. 

We conclude that the 7-parametrisation with a unique 
value of 7 can not be used for arbitrary companion 
masses. A small variation in the value of 7 covers the 
whole range of possible solutions — from unphysically 
expanded, to all possible physically acceptable configu- 
rations, and to mergers. As such, one cannot employ the 
7-formalism using 7 = 1.5 for an arbitrary binary, as 7 



Fig. 3. — Orbital angular momentum and energy for a 2Mq + 
1.5Mq binary (r d = 86 Rq and mj x = 0.38 Mq). Line styles are 
as in Fig.El 




lg J [G 1 / 2 M|/ 2 Ri/ 8 ] 



Fig. 4. — Orbital angular momentum and energy for a 2 Mq + 
0.35 Mq binary (r d = 86 Rq and m d ,x = 0.38 Mq). Line styles 
are as in Fig. [2] 

values must be fit for each system. To resolve the ap- 
parently negative «ce values required to explain the for- 
mation of many DWD systems, we propose that instead 
of a first common-envelope event, stab le mass transfer 
takes place. While it has been claimed (|Nelemans et al.l 
van der Sluvs et ahl 120061 ) that an initial phase of 
fully conservative MT is incompatible with the observed 
distribution of DWDs, this was demonstrated only for 
donors that started mass transfer near the beginning of 
the RGB. The degree of orbital expansion is greatly re- 
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duced if one introduces an inefficiency in the process — 
non-conservative MT. In this case, the second CE can 
proceed with physically reasonable values for acE- To 
avoid an initial CE, the first episode of mass transfer 
must proceed in a dynamically stable manner. 

4. STABILITY OF MASS TRANSFER 

Whether, once started, MT will proceed in a stable or 
unstable manner depends on the response of the donor 
and its Roche lobe to the MT. Mass loss is a perturbation 
that brings a star out of hydrostatic and thermal equi- 
librium. In order to re-establish these equilibria, the star 
will either expand or contract, first restoring hydrostatic 
equilibrium, and then, on a longer timcscalc, thermal 
equilibrium. Like the star's radius, its Roche-lobe radius 
also changes in response to the MT. 

As a stability condition, we thus demand that upon 
MT the new donor radius should remain within the 
donor's new Roche lobe. The dynamical stability of 
the MT depends on the adiabatic response of the stel- 
lar radius to mass loss. To consider stability in de- 
tail, we will use the li near st ability analysis f ol lowin g 
iHiellming fc Webbinkl (|1987fl ; ISoberman et al.l (jl997l ). 
We define the adiabatic mass-radius exponent (the 
donor's response to the mass loss on an adiabatic 
timescale) as 



ad 



d log r d 
dlog m d/ a 



(12) 



and the Roche lobe's mass-radius exponent (the Roche- 
lobe response to the MT) as 



RL 



d\og r R ,L 

G?l0g Tfld 



(13) 



The criterion for dynamical stability of MT then can be 
written as Cad > Crl- 

If this criterion is satisfied, the donor is able to re- 
cover its hydrodynamical equilibrium while remaining 
within its Roche lobe. It will then try to recover its 
thermal equilibrium on the Kelvin-Helmholtz timescale. 
The change in its equilibrium radius can be expressed 
with the use of the equilibrium mass-radius exponent 



c 



cq 



d log r d 
dlog irid 



(14) 



cq 



If, in addition to the dynamical-stability condition, the 
condition ( cq > Crl a l so holds, then MT will proceed 
on the nuclear-evolution timescale in a secularly stable 
manner. If Cad > Crl > Ccq, MT will be driven by ther- 
mal readjustment and will proceed stably on the thermal 
timescale (r t h). In the case of a red giant, the radius of 
the star in complete equilibrium is a function predomi- 
nantly of its core mass, and as such its thermal response 
is usually adopted to be Ccq ~ 0, though it can vary. We 
will now consider these responses in more detail. 

4.1. Adiabatic response 

The adiabatic mass-radius exponent for giants can be 
obtained cons idering the case of a conden sed polytrope 
with n = 3/2 (jHiellming fc Webbinld[T987l) . The value of 
Cad can then be found to within 1% accuracy using the 



formula from ISoberman et al.l (|1997l ): 



c 



2 M n 



ad ,, 



1 l -M r . 



3 1 - M c 3 1 + 2M< 



-0.03il/ c +0.2- 



l + (l-M c )-e' 
(15) 

where M c = m c (Mo)/m is the mass fraction of the he- 
lium core. 

4.2. Roche-lobe response 

We expand the logarithmic derivative of the orbital 
separation with respect to the donor mass as: 

(9 In a <91n(r RL /a) d In q 
Crl = H — — — — • (16) 



d In md 



d In q d\n to^ 



Here q — mr\/m a is the mass ratio. The first term, 
d In a/d In m, is solely due to the mass loss or MT and 
can be found using cq. ()lll) : 



d In a 
din md 



m d a 



(17) 



We will be most concerned with stability at the onset of 
a fast initial MT phase (see §4.3, §5.3). In this case we 
may assume that the MT rate from the donor exceeds the 
wind mass-loss rate significantly, so that for the purpose 
of this stability analysis we can assume that 5 = (Sec 
Then: 



d In a 
91nmd 



2m d 



2m a — md?7z a (l — (3) 



m a (m d + m a ) 



(18) 



where /3 is the fraction of the transferred mass that is ac- 
creted by the secondary star. We note that the presence 
of wind makes MT more stable, as the orbit always ex- 
pands in response to wind loss. As such, with 5 = 
we would obtain a stricter criterion for MT stability. 
The above equation demonstrates that the Roche-lobe 
response is a function of the MT conservation factor; in 
general, Crl ~ CrlC^q)- 

The second term in Eq.[T|)] consists of the Roche lobe's 
response to the change in mass ratio, which can be de - 
scribed using Eggleton's approximation (|Eggletonlll983|) : 



d\n(r RL /a) _ 2 _ q 1 / 3 l^q 1 / 3 + 1/(1 + q 1 ' 3 ) 
dlnq ~ 3 3 0.6g 2 /3 + l n (i + g i/3) 



(19) 



(see also ISoberman et al.lll997T) . and the response of the 
mass ratio to the change in donor mass 



d In q 
d In md 



1 



(20) 



We can then compute Crl from Eq. ll6l using the nec- 
essary parts from Eqs. lTSl IT91 andrSTil 

4.3. Stability boundary 

By comparing Cad an d Crl(/?) we can find a 
/3max(g, m d ,c), such that for all (3 < ^ max , Cad > Crl(/3) 
and therefore MT will be dynamically stable. We vi- 
sualize this in Fig. [5] for the case of a 1.2 M Q red giant, 
with unchanging mass during its evolution, and al.l Mq 
companion at the start of MT, demonstrating how this 
condition changes for increasingly evolved red giants. 
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Fig. 5. — Adiabatic and Roche-lobe mass-radius exponents for 
a red giant of 1.2 Mq with a companion of 1.1 Mq, as a function of 
the red-giant core mass. The primary mass is held constant during 
its evolution up to RLOF. Roche-lobe mass— radius exponents are 
shown for different cases of mass conservation. 
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Fig. 6. — Adiabatic and Roche-lobe mass-radius exponents for 
a red giant of 1.2 Mq with a companion of 1.1 Mq, as a function of 
the red-giant core mass. The primary has a metallicity of Z=0.03 
and is evolved with wind mass loss. Roche-lobe mass-radius expo- 
nents are shown for different cases of mass conservation. 

It can be seen from Fig. [5] that in all the cases of 
fully non-conservative MT (fi = 0), the MT is dy- 
namically stable. We also note that fully conserva- 
tive MT is dynamically stable for cores with masses 
m d,c ^ 0.63 Mq, and in this case MT will proceed on 
the thermal timescale. 

In Fig.[S]we see that in the more realistic case in which 
we account for mass loss (including wind loss), as well 
as z = 0.03, our range of stability is extended even fur- 



ther, to the point of allowing fully conservative MT for 
large core masses (M<j iC > 0.46 Mq). With decreasing 
mass RLOF-drivcn MT is increasingly stable, e.g. for a 

1.1 + 1.0 M Q RG-MS system conservative MT is stable 
for M c > 0.42 M . 

Of primary interest to us at the present moment are red 
giants that are potential progenitors of the inferred older 
companions of observed DWD systems, hence our focus 
on systems in which the donor core mass reaches ~ 0.35— 
0.46 M (before the He flash). Note that RLOF-driven 
MT may often run on the nuclear timescale, such that the 
core mass can grow significantly during this phase. Hence 
we must consider those systems in which MT begins with 
a core mass as low as 0.25 M Q and up (see §5.3.2). 

Using the above condition for stable MT, we can define 
an upper boundary on the mass ratio for any given core 
mass of the donor in order for the initial MT phase to 
proceed in a stable manner. This upper boundary varies 
with conservation factor through Crl's dependence on (3 
(see Eq. [l6|) . For the range of core masses of interest here, 
we find that if MT occurs with a conservation factor no 
greater than /3 ~ 0.3 — 0.5 (depending on initial mass 
ratio), it will be dynamically stable. 

We also note that the stability criterion we use may 
not be final, as it is based on values of Cad for condensed 
polytropes, not for real giants. As far as we know this is 
a more restrictive criterion than the ones typically used 
in the literature, which are frequently based on the mass 
ratio. For example, f or giant donors q = 1.2 i s considered 
to be a threshold in iBelczvnski et all (|2008l see the ref- 
erences therein). This condition is derived from detailed 
MT calculations resulting in runaway rates. Adopting 
the criterion based on the comparison of Cad ^ Crl(/3), 
we are at the most conservative limit, and even this limit 
predicts that some systems can have stable, fully conser- 
vative MT. We also check when the MT rate grows ex- 
ponentially (indicating instability) in our MT sequences 
and compare it to the predicted value of /3 ma x- 

4.4. Thermal equilibrium response and the end of 
thermal-timescale mass transfer 

In the case of a red giant, the radius of the star in com- 
plete equilibrium is a function predominantly of its core 
mass, and its thermal response is usually considered to 
be Ceq ~ 0. However, the detailed comparison of giants 
of different total masses, but with identical core masses, 
indicates a dependence on the tota l mass as well (see e.g. 
Fig. 1 in Ivan der Sluvs et al.ll2006h . In particular, for gi- 
ants with masses < 1.2 M Q , £ eq is negative and a function 
of the giant's core mass. To find it more precisely, one can 
build a sequence of stellar models in thermal equilibrium 
with the same chemical composition and the core mass, 
starting with 1.2 M Q and then decreasing the mass. This 
is done by imposing a very slow mass loss on a 1.2 Mq 
giant and switching off the chemical evolution between 
the relaxed stellar models. 

We find that a giant with a 0.37 M Q core has Ccq ~ 
—0.4 as long as the giant mass is between 0.6 M Q and 

1.2 Mq. For masses below 0.6 Mq, £ e q becomes positive. 
We know that as long as £rl > Ccq, MT proceeds on the 
thermal timescale, and it will switch to nuclear-timescale 
MT when Crl < Ccq- We can then find the mass ratio 
g cl -it when the condition for the thermal-timescale MT 
(TTMT) is no longer satisfied. For instance, we find that 
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if a 1.2 Mq giant with a 0.37 Mq core has a companion of 
1.1 M Q and MT proceeds with /3 = 0.3, then Crl = Ccq 
when the mass ratio in the system becomes 0.75 (i.e., 
when the donor mass is decreased to ~ 0.9 Mq). Due 
to some inertia in the star's response, the TTMT may 
proceed afterwards, but no longer than for ~ tth, and at 
g crit sa 0.75 the MT rate is likely to be at its maximum. 
In the case of a 1 Mq companion, the maximum TTMT 
rate is expected to occur at g cr it ~ 0.85 (donor mass 
is decreased to 1M©). Though it does not yet give us 
strong predictive power on when the TTMT must stop, 
it helps us to understand the results of our simulations 
qualitatively. 

5. THE FIRST PHASE OF MASS LOSS 
5.1. Stellar- evolution code 

In order to evaluate our model we perform numeri- 
cal calculations using the binary stellar-evolution 
code originally developed by E ggleton ([Eggleton I1971L 
119721: lYakut fc Eggleton! 120051 and referen c es th erein) 
and updated as des cribed in IPols et al~l (119951 ) and 
iGlebbeek et al.1 (|2008l ). The code solves the equations 
of stellar structure and evolution for both components 
of a binary simultaneously using an implicit scheme over 
an adaptive mesh. We model each star using a grid size 
of 200 mesh points, as this consistently provides sta- 
ble results. This results in the model being quite sta- 
ble against very short-timescale instabilities, allowing us 
to quickly evolve our models up the asymptotic giant 
branch (AGB). Simultaneous calculation of both compo- 
nents is essential in order to account for non-conservative 
mass transfer, as ea ch star's evolution is no longer separa- 
ble from the other (lEggleton fc Kiseleva-Eggletonll2002l : 
lYakut fc Eggletonll2005D .~ 

Op acity tables are taken from OPAL (Iglesias et al. 
199a. with the low- t emper ature range taken from 
Alexander fc Ferguson! ([1994ft . A diffusion equation 



models convective mixing for each of the composi- 
tion variables, and overshooting is modelled as in 
iSchroder et al.l (jl997j ). The helium flash in the degen- 
erate core of low-mass stars is avoided by substituting 
the stellar model with one in which helium has just been 
ignited in the core. The initial metallicity for all models 
is assumed to be roughly solar (X = 0.70, Y = 0.28, Z = 
0.02). 

We account for both RLOF and a stellar wind in the 
bou ndary conditions for t he mass, assuming a Reimers- 
like CE eimersl 119751 11981D model for stellar- wind-driven 
mass loss: 



M = 4.0 x 10 -13 



n 



R L Mq _ x 

M(7)yr , 

RqL q m oy 



(21) 



where we set 77 = 0.2. 

Egglcton's code provides two means of parametrizing 
the mass loss from the donor during RLOF. One way is 
to allow the donor's radius to slightly exceed its Roche 
lobe, whereupon the mass loss is computed in terms of 
this excess: 



CMS x 



log 

mL 



(22) 
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2 The current version of ev is obtainable on request from 
eggletonl011nl.gov, along with data files and a user manual. 



Fig. 7. — Mass-transfer rates according to different prescriptions, 
for a 1.2 + 1.1 Mq RG-MS binary, with a period at the start of 
RLOF of 100 d. The solid and dotted lines indicate the CMS for- 
mulation, the dashed line is for the CMT formulation. A mesh of 
200 points is used in each case. 



Here CMS is an arbitrary parameter; a larger value is ex- 
pected to provide more appropriate mass-transfer rates, 
as this allows for a smaller Roche-lobe overflow through- 
out the mass-transfer evolution, and hence smaller time 
steps taken by the code. A lower value may provide MT 
rates that are too small for the case of TTMT, though the 
donor is exceeding its Roche lobe more than in the case 
of a bigger value, and the code takes larger time steps. 
In Fig. [3 we plot the evolution of the mass-transfer rate 
for a 1.2 + 1.1 Mq binary with an initial period of 100 
days while varying CMS. We find that so long as a rea- 
sonable value is chosen (i.e. one for which runaway MT 
is not induced artificially), the same salient evolutionary 
features are seen, albeit at slightly different moments in 
the MT history. Notably, in all cases MT enters a fi- 
nal phase of slow, steady nuclear-timescale MT following 
a thermal-timescale phase and a pause (see below). In 
addition, for each case final values of m^^ and P agree 
within ~ 1%. We adopt CMS = 10 as our default value 
— it is both found to be the most appropriate value in 
our calculations, and is generally recommended for donor 
stars in the mass range considered here (Eggleton, pri- 
vate communication). 

An alternative is a brute-force formulation, labelled in 
the code as CMT, in which the potential is calculated 
at each grid point within the donor star. This is then 
used to calculate mass flux £ away from the donor as 
a function of depth, zero below the L\ surface (again, 
the donor radius is allowed to slightly exceed its Roche 
lobe). The MT rate is then calculated as M = CMT • £, 
where CMT is an arbitrary parameter. This method is 
designed for contact binaries. In Fig. [71 we plot the MT 
evolution for the same system using this prescription: 
for CMT = 1, we find a MT rate history comparable 
with that obtained for CMS = 10, but with more rapid 
oscillations. At higher values we find runaway MT which 
we again assume to be artificial, as they seem to be the 
result of the arbitrary coefficient . 

The oscillatory pulsations noted above, when a result 
of numerical error, arise as an artefact of defining too 
large mass shells when implementing the mesh used to 
compute our model, partly due to our choice of a rela- 
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tively small number of mesh points (200). This causes 
the boundary between the star's convective envelope and 
the core to vary rapidly at high MT rates, resulting 
from overly large mass shells changing back and forth 
between a convective and a non-convective state at the 
core boundary. These pulsations appear to be much 
stronger when using the CMT prescription (Fig. [7]), sug- 
gesting greater instability. However, these pulsations do 
not have a strong influence on the evolution of our model 
and its outcome, and for our purposes may be ignored. 

5.2. A maximum conservation factor 

Whether the dynamical-stability condition Cad > Crl 
is satisfied or not depends on (3: Crl is dependent on 
/3 through the responses of both the mass ratio and the 
orbital separation to MT (see Eqs.[IS[l8]and[l9]). As Cad 
and Crl are both monotonic functions of mass loss, it is 
a sufficient condition that MT be stable at the onset of 
RLOF in order to ensure stability throughout. We reject 
system configurations in which this condition is violated. 

We realise that there is a difference between a true 
dynamical stability (or instability) — the one arises on 
time-scales shorter than dynamical-timescale, before a 
star obtain hydrostatic equilibrium — and the dynam- 
ical stability used commonly for the purpose of studies 
of mass transfer stability. The latter one, as described 
earlier in §4, operates on timescales longer than the dy- 
namical timescale (as all the analyzed stellar model arc 
in hydrodynamical equilibrium) but much shorter than 
thermal timescale, so the entropy of the stellar layers is 
not changed. The code we use always generating stel- 
lar models in hydrostatic equilibrium, and as such is not 
capable of modelling a true dynamical-timescale insta- 
bility, however is able work on a timescale much shorter 
than is required to change the entropy. As such, through 
out mass transfer sequences, when the rate of the MT 
exceeds TTMT, we efficiently obtain Cad as discussed 
above, but for actual stellar models (instead of composite 
polytropcs). We find that with (3 close to the maximum 
stable value, we obtain MT rates exceeding the thermal- 
timcscale MT rate, although this is achieved only after 
~ r t h has passed after the onset of RLOF. 

As an example, we model the case of a 1.2 M© red gi- 
ant with a 1.1 M companion, with an initial period of 
100 days, which initiates MT at a core mass of approx- 
imately 0.345 Mq and with a donor radius of 46.4 i? . 
For a conservation factor greater than /3 max ~ 0.89, run- 
away MT ensues (see Fig. [8]): MT stability is unable to 
recover on the thermal timescale for ft — 0.9, leading to 
the code crashing. Wc interpret this as indicative of a 
runaway MT event. In this case the MT rate approaches 
and finally exceeds the thcrmal-timcscalc-MT rate. Note 
that the maximum conservation factor found here is con- 
siderably higher than that found when taking the donor 
as a condensed polytrope considered with the core and 
the total masses as in the described above giant, in which 
case /3 

max 

ps 0.32. 

5.3. Detailed evolution 

Note from Fig. [7] that the MT phase appears to pass 
through two distinct sub-stages. Using the example 
above (a 1.2+1.1 Mg progenitor system with Pi = lOOd) 
for the case (3 = 0.3, we observe that there are in fact 
three such stages. 
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Fig. 8. — Mass-transfer rates for different conservation factors /3 
in a binary of ; = 1.2 Mq, m a ; = 1.1 Mq, and with an orbital 
period at the beginning of RLOF of 100 d. 
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Fig. 9. — Time evolution of the MT rate (dashed line) and mass 
ratio (solid line) after the onset of MT. The MT proceeds on the 
thermal timescale, and Qcrit ~ 0.75. 

At the onset, the MT rapidly accelerates to a timescale 
comparable to the donor's thermal timescale, reaching 
a dramatic peak M before rapidly falling off. This 
turn-off occurs approximately r t h after the condition 
Crl(<?) < Ccq is met (see Fig. [9]), as the mass ratio reaches 
<7crit- This initial stage strips the envelope down (remov- 
ing ~ 50 — 70% of the envelope mass), to which the donor 
reacts by contracting until the binary becomes detached 
(underfilling the Roche lobe by ~ 5%). This results in a 
pause in MT lasting for ~ 1 — 4 Myr (this recovery time 
is approximately the thermal time of the whole donor, 
whereas the thermal time of only its outers layers played 
a role in the determination of the TTMT rate), during 
which time the donor's core grows slightly until RLOF 
resumes. 

This pause is followed by another stage of MT on the 
star's nuclear timescale, in which the remainder of the en- 
velope is transfe rred, save f or a final ~ 10~ 3 M . While 
often neglected (jHanlll99"8f l. in the course of the pause 
and nuclear-timescale phase of MT the donor's core mass 
can increase by ~ 20 — 30%, having a profound effect on 
the outcome of the system. 

As for the donor's remaining envelope mass, at the end 
of MT this collapses (on its thermal timescale) onto the 
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Fig. 10. — Evolution of the donor radius (solid line) and MT 
rate (dotted line) as a function of donor mass during MT, for an 
initially 1.2+1.1 Mq binary with an initial period of 100 days. The 
Roche-lobe radius is denoted by dashed blue line. 



Fig. 11. — Period growth with core-mass evolution for a 1.2 + 
l.lMg binary with an initial period of 100 d, broken down in the 
three stages of the first MT phase. 



surface of the degenerate core (iDeinzer k, von Sengbuschl 
[T^IIvanovdl20nt[jvrsthainll2011l) . The same behavior 
is observed in our model, and the amount that remains 
in the envelope at the end of MT is generally ~ 3% of the 
original envelope mass. This means that for smaller core 
masses, we see thicker envelopes remaining. However, 
using our code it is difficult to make a precise estimate 
of the final mass of the envelope remnant, as our model 
breaks down at the collapse of the envelope for low core 
masses. For the purpose of this study, we assume that the 
MT ends very shortly after the code breaks down, and 
that the remaining mass in the envelope will be burned, 
i.e., we take the total mass at the end of the MT phase 
as our final WD mass. We note that this may introduce 
a slight bias in the results, pushing the final mass ratio 
somewhat closer to unity and the final period to very 
slightly shorter values in the least massive systems. 

In order to better illustrate each stage of the first MT 
phase, we display the evolution of the MT rate and donor 
radius as a function of the donor mass in Fig. llOl and 
the evolution of the orbital period as a function of the 
donor's core mass in Fig.QTJ Figure IT2"1 shows the donor's 
radius as a function of its core mass, and we note that 
the radius of the donor during the nuclear-timescale MT 
is significantly larger than what we would expect of an 
unperturbed star. 

To investigate how the above evolution varies under 
different initial parameters, we now recompute the evo- 
lution of our 1.2 + 1.1 Mq binary with P; = 100 d, and 
vary the conservation factor (Fig. [T3"|) . initial donor mass 
(Fig. Q3]), and initial mass ratio (Fig. [T5|) . This exercise 
demonstrates that for initial conditions leading to the 
first phase of MT, 

• the conservation factor is anti-correlated with the 
final period, 

• a larger initial donor mass (for constant initial mass 
ratio) leads to a longer final period, 

• the initial mass ratio is anti-correlated with the fi- 
nal period. 

In the next section we will also vary the initial pe- 
riod, which will naturally lead to larger core masses for 



cr 
cr 




0.15 0.2 0.25 0.3 0.35 0.4 0.45 



Fig. 12. — Deviation from the core mass-radius relation for the 
donor star from the system in Fig. llOl Solid line: unperturbed star, 
dashed line: perturbed donor star. 
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Fig. 13. — Evolutionary tracks for a 1.2 + 1.1 Mq binary with 
Pi = 100 d through stable, non-conservative MT, for fi = (top), 
0.2, 0.4, 0.6, and 0.8 (bottom). 
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Fig. 14. — Evolutionary tracks for an initial binary system with 
P = 100 d through stable, non-conservative MT, for initial = 
1.15 Mq (bottom), 1.20 Mq, 1.25 Mq (top), and a constant mass 
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Fig. 15. — Evolutionary tracks for an initial binary system with 
Pi = 100 d through stable, non-conservative MT, for q\ = 1.200 
(bottom), 1.143, 1.091, 1.043, and 1.008 (top). 

larger initial separations. In addition, we will compute 
the further evolution of such binaries through the RLOF 
of the (former) secondary, which — due to the dramatic 
change in mass ratio — will proceed in an unstable man- 
ner, leading to a typical CE. 

6. THE SECOND PHASE OF MASS LOSS 

6.1. An Ensuing Common Envelope 

In the intervening period following the end of the first 
phase of MT, further orbital expansion may be driven by 
a wind from the (former) secondary (i.e., the originally 
less massive component). However, as noted above, in 
many cases our code has difficulty modelling the collapse 
of the envelope at the end of stable MT, and therefore 
cannot progress through this intermediate stage. Hence, 
in continuing to the second phase of mass loss we ignore 
wind loss from the secondary. Those models which do 
survive the end of the initial MT phase suggest the rel- 
ative error introduced in the intermediate period (P m ) 
by ignoring the wind-driven expansion is on the order 
of 5 — 10%. For the subsequent CE phase, the resulting 
error in the radius at RLOF is small, as at the onset 
of RLOF the radius will only be related to the period 
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Fig. 16. — Full evolutionary tracks for a binary with initial masses 
of 1.2 + 1.1 Mq through an initial stage of stable MT followed by a 
CE, for Pi = 1, 5, and 25 d with respect to donor mass. Circles indi- 
cate remaining (core) mass of primary and orbital period after the 
system has gone through the initial phase of stable MT. Arrow end- 
points indicate the final (WD) mass of the secondary and orbital 
period after the second, unstable episode of mass loss (CE phase). 
Hence, an arrow slanted to the right indicates a final system where 
the secondary WD is more massive than the primary. Black lines 
trace the evolution of the system through the first (dashed) and 
second (solid) phases of mass transfer. Black stars (from doublc- 
lincd binaries, table 1) and triangles (from Kilic ct al. (2010)) mark 
parameters of observed DWDs (orbital period and mass of the in- 
ferred secondary, see Tables [l] [2]| . 



weakly through P 2 ~ a 3 ~ R 3 . Since inverting the m c - 
R relation gives the core mass a roughly logarithmic de- 
pendence on the radius, we may take the error in the 
final core mass of the secondary at the onset of the CE 
as being negligible. As the post-CE separation is directly 
proportional to the initial separation in the q;ce formal- 
ism, we may take the relative error in Pp os t-CE as being 
roughly the same as that in P m . 

From this point we can easily estimate the resulting 
final binary configuration for a given system after a sec- 
ond episode of mass loss in which a common-envelope 
phase ensues (assuming that the acE formalism holds, 
and computing the orbital change directly from the bind- 
ing energy in our models). We use a core- mass ra- 
dius relation from the parametrized stellar models of 
(jHurlev et al.l [2000) in order to estimate the core mass 
of the secondary at which the onset of a CE takes place, 
as this approximates the radii of all our models to within 
a few percent. As a CE event is a dynamical-timescale 
process, we neglect core-mass evolution over the course 
of the CE phase. 

As an example, we evolve a set of 1.2 + 1.1 Mq systems 
with /? = 0.3 and varying initial period as above through 
a second, unstable episode of MT, plotted in Figs. [Hal 
1171 These form a fairly linear track of end products in 
both m — log P space and q — log P space, with the fi- 
nal periods, mass ratios, and core masses of the inferred 
younger companions matching nicely with observations. 
Initial period aside, the final period is most strongly de- 
pendent on the efficiency of the CE, the lower values of 
acE driving the final separation to smaller values. Fur- 
ther modulation of these results may come from varying 
the initial donor mass, initial mass ratio, and conserva- 
tion factor. 
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Fig. 17. — Full evolutionary tracks for a binary with initial 
masses of 1.2 + 1.1 Mq through an initial stage of stable MT fol- 
lowed by a CE, for Pi = 1, 5, an d 25 d, plotted against mass ratio. 
Lines and symbols are as in Fig. 1161 Due to the much greater un- 
certainty in the masses of the companion WDs in Table [2] we plot 
only those observed systems listed in Table [T] 



6.2. An Exception: WD 1101+364 

Among those double-lined DWDs in which neither 
component's progenitor has undergone the He flash, we 
determine then that the stable MT+CE channel pro- 
vides a reasonable explanation for the formation of those 
systems with q( inferred older/younger) < 1. This ex- 
cludes only one system, WD 1101+364 (Fig.[T71) whose 
mass ratio is greater than unity, implying the orbit 
contracted with the loss of the initial donor's enve- 
lope. Attemp ts to model this system as the product of 
a dou ble- CE (jNelemans et alJl2000t Ivan der Sluvs et al.l 
20061 ) found this could only be done with an unphysi- 
cal value for aX for the first mass loss phase; this was 
part of the motivation for t he de velopment of the 7- 
formalism. iNelemans et al.l (|2000D suggested this sys- 
tem could better be modelled as the product of two 7- 
prescribed events, with the initial mass inferred from this 
reconstruction implyi ng a differenc e in c ooling ages of 
- 800 Myr. However. iMaxted et al.l (|2002t) found an ob- 
served difference of ~ 215 Myr, certainly irreconcilable 
with the predicted value. As they are quick to point out 
though, cooling models for helium white dwarfs remain 
quite uncertain, and drawing strong conclusions regard- 
ing their relative ages remains a difficult task. 

Compounding the problem of determining WD ages 
is the prospect of rejuvenation, as even the accretion 
of interstellar matter may pose a significant problem 
l|Nomoto fc Sugimotolll977f ). It is therefore quite pos- 
sible that even our understanding of which component 
of any DWD is older may be incorrect. Indeed, WD ages 
are based on models of cooling tracks whose uncertainties 
are highly dependent on the size of the remaining hydro- 
gen envelope, with thicker envelopes working as 'blan- 
kets' and causing WDs to appear hotter, hence younger. 
WD 1101+364 has the smallest component of any of the 
companion masses, and as we have seen above (§5) this 
corresponds to a thicker anticipated envelope. Such a 
turnaround in the respective age difference is not impos- 
sible. We note that a reversal of the mass ratio would 
place WD 1101+364 almos t directly o n the line traced 
by our example track, and iHanl (|1998f ) models this sys- 



tem as originating through a stable MT+CE scenario 
under similar assumptions, though this is hardly defini- 
tive. As well, we have neglected the effect of hydrogen- 
shell flashe s after the collaps e of the envelope here (e.g. 
Sect. 5.3 in lYoon et aLl[200l . which may reduce the en- 
velope mass significantly due to nuclear burning, mass 
transfer and mass loss. Shell flashes occur earlier and 
are stronger for lower-mass proto-WDs. 

Further model ling of WD 1101+364 in 

Ivan der Sluvs et al.l ( 20061 ) produced better results 
by varying the means by which the binary shed angular 
momentum. Assuming that in the first phase mass is 
lost directly from the donor in an isotropic wind, while in 
the second phase it is re-emitted from the accretor, gave 
a small set of solutions which were nearly within error. 
Yet a physical explanation for such processes to occur 
on a timescale short compared to the nuclear-evolution 
timescale, as assumed^ remains as of yet elusive. At 
best, we can say that something appears to be amiss in 
our understanding of WD 1101+364. 

7. DISCUSSION AND CONCLUSIONS 

In this paper we addressed the problem of doublc- 
whitc-dwarf formation, beginning with an analysis of the 
applicability of the previously proposed 7-formalism for 
an arbitrary binary. Originally, the 7-formalism was de- 
veloped in order to describe a dynamical-timescale phase 
of mass loss in terms of a single-valued parameter for 
all possible binaries, without a specific assumption of 
the physical characteristics of the evolution. This value 
was then found based on the study of currently observed 
DWD systems. 

We have demonstrated that, when parametrizing the 
consequent angular momentum loss with a single value 
(7), two dynamical-timescale phases of mass loss can only 
model the formation of arbitrary DWD binaries with the 
use of different values of 7 which arc very finely tuned for 
every mass ratio. The underlying physics of a first phase 
of dynamical-timescale mass loss in such a process, in 
which the envelope of the initial primary must be ex- 
pelled without significant shrinkage (or even widening) 
of the orbit, remains without an explanation under the 
formalism of a 7-prescribed envelope removal. However, 
the apparent violation of energy conservation during the 
first phase of mass loss as described by the 7-formalism in 
many cases indicates that the envelope ejection cannot 
be driven by any dynamical process in such instances, 
but instead appears to be mimicking the stable, non- 
conservative MT that runs on a longer timescale — ther- 
mal or nuclear. 

Indeed, in §5,6 we demonstrate that such stable, non- 
conservative MT provides a physically motivated, easily 
implemented mechanism in lieu of any paramctrization 
of the energetics or angular momentum balance during 
the first mass loss phase, which here obscures the phys- 
ical interpretation and, as we have seen, requires care- 
ful fine tuning. Circumventing the need for any singlc- 

3 Note that while van der Sluvs et al. (2006) claim to assume 
that envelope ejection described by the 7-prescription must occur 
on the dynamical timescale, they in fact make a looser assumption; 
for the mechanism to work, ejection of the envelope on a timescale 
shorter than the nuclear-evolution timescale — so that the core 
mass doesn't change — is sufficient, and this is their actual as- 
sumption. 
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valued paramctrization, wc compute the evolution ol any 
arbitrary low-mass DWD progenitor for which the in- 
ferred older component has not undergone the helium 
flash through the self-consistent computation of the evo- 
lution of a binary through stable MT (assuming the sta- 
bility conditions outlined in §4, 5 hold), followed then by 
a CE. 

Hence, this scenario allows us to reproduce the period 
and mass ratio of the observed DWDs with low-mass 
companions (M < 0.46 M©) with only a single combina- 
tion of initial binary-component masses (1.2 + 1.1 M©). 
This is very promising, but more work is needed to deter- 
mine whether we can indeed explain all the observed low- 
mass DWDs in which the older component is less mas- 
sive, including the differences in cooling age using this 
scenario. The latter, of course, would require stronger 
constraints on WD cooling tracks (including the effect of 
rejuvenation) than exists now. 

For MT evolution, we analyze how non-conservative 
MT should be in order to proceed stably. We did not 
discuss in detail the exact mechanism for expelling mass 
from the system during MT, and do not propose what the 
exact value of the efficiency should be, but rather study 
the consequences of non-conservation during MT with 
an adopted conservation factor. We may argue however 
that MT is rather likely to be non-conservative, espe- 
cially during the TTMT phase. The strongest argument 
for this is due to the very high MT rate expected: some 
of the binary systems, at the peak of TTMT, could have 
MT rates above the Eddington limit (~ 10 3 M Q yr _1 
for MS accretors in the mass range we considered). Ob- 
servationally, though, it is hard to make any constraint: 
the TTMT phase is of course very short lived. We note 
however that among observed binaries undergoing MT 
with giant donors, e.g., in symbiotic stars, there are ob- 
served systems showing powerful jets — clearly, the MT 
is not fully con servative the re (e.g., in sy mbiotic stars 
like CI Cygni (lliiimal U982T) . CH Cygni dTavlor et al.1 
EM)), even though the MT proceeds there on a nu- 
clear timescale. Young pre-MS stars are the only ob- 
served non-compact objects that accrete at rates almost 
comparable to our TTM T rates, and are kno wn to have 
jets as well e.g., HH 30 (|Anglada et al.ll2007t) . Thus our 
assumption that bipolar re-emission from the accretor 
drives non-conservation in MT, and therefore the spe- 
cific angular momentum of the lost material is that of 
the accretor, seems at least in principle quite reasonable. 

In the context of earlier studies, it is important to 
note that the progenitor donors with M > 2 M© of- 
ten considered previously would not be able to form a 



WD immediately following an episode of MT — they 
would instead form a low- mass He star with a lifetime of 
up to 0.5 Gyr (|Yung clson 20081). relevant to the age dif- 
ference of many of t he so lutions described in table 5 of 
Ivan der Sluvs et al.1 (j2006[ K In this study we find that in 
the low-mass regime the stable MT + CE channel, with 
an initially 1.0—1.3 M© donor, provides a natural expla- 
nation for all but one of the observed double- white-dwarf 
binaries for which the inferred older component is below 
0.46 M Q . 

Thus, systems such as PG1115+116, in which the 
younger companion is the more massive, can easily form 
as a result of the orbital expansion during the initial 
phase of RLOF, followed by a CE. The relatively small 
delay between the two phases of mass loss from the sys- 
tem, inferred from observations, is also easily justified 
by the accelerated evolution of the accretor as a result of 
the first MT episode. Of course, a strong determination 
of the relative ages in a DWD is a difficult task, further 
confounded by the possibility of rejuvenation. 

The system WD 1101+364 presents a clear challenge 
in that its older component is the more massive, yet it 
cannot be explained as the product of a double com- 
mon envelope without resorting to unphysical efficien- 
cies in the first phase of mass loss. One possibility may 
be some intermediary pr ocess, such as that suggested 
by iNelemans et al.l (|2000ft , in which the envelope of the 
initi al primary is lost wi t hout significant orbital shrink- 
age; Ivan der Sluvs et al.l (|2006| ) suggests this stage may 
be understood as a phase of very rapid wind loss. How- 
ever, an obvious mechanism for this remains lacking. It 
is also possible that the formation of a small circumbi- 
nary disk in the first phase of mass transfer may allow 
for some shrinkage of the orbit, however there is a lack 
of evidence at the present moment to support such a 
hypothesis. Perhaps the most likely explanation, given 
the small secondary mass in WD 1101+364 and there- 
fore its thicker hydrogen envelope, is that the resulting 
uncertainty in the cooling tracks means that the 0.29 M© 
companion is older than previously thought. 
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TABLE 1 

Double White Dwarfs with Known Masses and Periods 





System 


Ml[M ] 


M 2 [M Q ] 


Period [d] 


q 


1 


WD 0136+768 


0.35 ± 0.02 


0.46 ± 0.03 


1.407227 


0.793 


2 


WD 0957-666 


0.32 ± 0.02 


0.37 ± 0.02 


0.06099 


0.885 


3 


WD 1349+144 


0.35 


0.44 


2.2094 


0.793 


3 


WD 1349+144* 


0.44 


0.44 


2.2094 


0.793* 


4 


WD 1101+364 


0.35 


0.29 


0.14458 


1.149 



Note. — Known double-li ned DWD systems with established orbit al pa rameters, for which the (inferred) olde r WD (Mi) is less than 
0.46 M p,. Results taken from IBragaglia et al.l 1(1990 ) [2] . IMarshl l[T995l ) [4] . IMoran eTaLl lfl997l ) [2]. IMaxted et al.l (f200^) [1.2]. IKarl eFaTI 
(2003) [3]. For WD 1349+144, we note that the measured mass of Mi is 0.44 M (marked with an asterisk), however based on the (much 
more trusted) observed q- value of q a 0.793, Mi ?s 0.35. We plot the latter. Also note that here we define q = (inferred older WD)/(infcrrcd 
younger WD), opposite to the usual convention. 



TABLE 2 

Double White Dwarfs with Partially Known Masses and Periods 



System 


A/i, min [M ] 


Mi, max [M Q ] 


M 2 [M ] 


Period [d] 


J0022+0031 


0.35 


? 


0.38 


0.492 


J0022-1014 


0.19 


0.45 


0.33 


0.079 


J 1234-0228 


0.09 


0.22 


0.23 


0.092 


J1625+3632 


0.08 


0.17 


0.20 


0.233 



Note. — Known single-lined DWD systems with established period and mass of one component. Where possibl e, upper an d lowe r 
bounds on the companion masses are taken from the source paper (to varying confidence levels). Results taken from iKilic et all H2010I ). 
The observed white dwarf is presumed to be the younger one (M2). 



